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adiabatic passage. These opportunities are associated with pulse excitation of ladder- 
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off-resonant infrared radiation sweeps the laser-induced Stark-shift of a two-photon 
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duced by the first ultraviolet laser. A judiciously delayed third pulse scatters at this 
coherence and generates short-wavelength radiation. A theoretical analysis of these 
problems based on the density matrix is performed. A numerical model is developed 
to carry out simulations of a typical experiment. The results illustrate a behav- 
ior of populations, coherence and generated radiation along the medium as well as 
opportunities of efficient generation of deep ( vacuum ) ultraviolet radiation. 
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1 Introduction 

A technique of nonlinear-optical frequency conversion in gaseous media provide a 
well-established powerful tool for the generation of tunable coherent short-wavelength 
vacuum-ultraviolet (VUV) radiation [1]. Many applications have been found for 
these light sources especially in spectroscopy. However, due to the relatively small 
off-resonance non-linear susceptibilities and severe problems associated with absorp- 
tion and dispersion in resonance atomic or molecular gases conversion efficiencies 
obtained with this technique is rather small, typically in the range of 10~ 6 -10~ 4 . 
(For a review, see e.g. [2, 3]). Negative accompanying resonance effects can be 
reduced by manipulating the coherently-driven medium through nonlinear inter- 
ference effects at quantum transitions (see e.g. the review in [4, 5, 6]) in order to 
minimize resonant absorption and maximize non-linear susceptibilities. One of the 
manifestations of destructive interference at quantum transitions is electromagnet- 
ically induced transparency (EIT) [7, 8] and coherent population trapping [9], which 
lead to the cancellation of absorption for some of the coupled fields and may im- 
prove phase matching. By the appropriate choice of detuning and intensities, VUV 
radiation that is generated in a frequency conversion process has been produced 
with conversion efficiencies up to 10 -2 [10, 11, 12]. In strongly-driven media, one 
cannot think about nonlinear-optical response in terms of susceptibilities in their 
original sense. While EIT may lead to increased conversion efficiencies, the use of 
large coupling fields to produce transparency in a Doppler-broadened medium is in 
conflict with the requirements for maximizing the nonlinear susceptibility. Despite 
the constructive interference, the large Autler-Towncs splitting results in much re- 
duced nonlinear susceptibilities compared to resonant values. Furthermore, when 
tunability is required, resonances with the generated field cannot be used to provide 
an enhancement in the nonlinear optical response of the medium. 

In a two-level atom, optical polarization depends on the coherence between the 
ground and excited states, and reaches a maximum of 1/2 for an equal amplitude 
coherent superposition of ground and excited states. Therefore, in a medium com- 
prised of a high density of dark-state atoms, a relatively-large nonlinear polarization 
can be produced while absorption for the driving resonant fields is decreased. The 
system may then be regarded as an oscillator driven with maximum amplitude. If 
another electromagnetic field is introduced, the beat frequency produced by inter- 
action with the oscillator may be generated with high conversion efficiency, which in 
principle can approach unity. In this system the nonlinear polarization responsible 
for generation may become as large as the terms responsible for absorption. Thus, 
nonlinear interference processes in media with maximum coherence can substantially 
improve frequency conversion between the coupled fields. Frequency-mixing pro- 
cesses at maximum coherence have been studied theoretically [13, 14, 17, 18, 15, 16] 
as well as demonstrated experimentally [13, 14, 19, 20] both with adiabatic laser 
pulses and for steady-state conditions. In the latter case, the lowest relaxation rate 
and therefore maximum coherence is achievable at the Raman transitions associated 
with the ground electronic state. Consequently, in four-wave mixing experiments, 
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maximum coherence has been established in a A- type coupling scheme, similar to 
that employed for coherent anti-Stokes Raman scattering (CARS). A-type coupling 
gives limits to the shift of the generated radiation to the short-wavelength ranges 
compared to the ladder-type schemes. Besides that, laser-induced Stark shifts, 
which are intrinsic to strong coupling, may perturb the adiabatic population dy- 
namics and prohibit the preparation of maximum coherence [22, 23, 24]. 

When the pulse duration is much shorter than the coherence relaxation rates 
in a quantum system, coherent quantum control of populations based on Stark- 
chirped rapid adiabatic passage (SCRAP) can be performed that does not require 
Raman type ground-state coherence. Physical principles of this promising technique 
and proof-of-principle experiments are described in a recent publication [25]. In 
SCRAP, two-photon (in the general case, multi-photon) detuning of the driving 
field is controlled through a Stark shift of the upper level by an auxiliary strong 
laser pulse Est(^st), which is time-shifted with respect to the pump laser. The 
process is closely related to rapid adiabatic passage controlled by a frequency chirp 
in the pump laser [26, 27]. 

This chapter investigates potentials and features of deep (vacuum ultraviolet) 
generation at maximum coherence based on the SCRAP technique [25]. However, 
instead of maximum population transfer from the ground to an excited state under 
the driving field E\(w\), the emphasis is placed on generating persistent maximum 
coherence at the frequency 2wi. This is related to robust transfer of half of the 
population from the ground to the two-photon-resonant upper state. This process 
is substantially effected by a time-dependent Stark shift of the two-photon resonance 
induced by both control (auxiliary) and fundamental (pump) time-shifted pulses. 
A third probe field {0J2) scatters at oscillations at 2a>i and generates difference- 
and sum- frequency deep UV radiation at lo t = 2cji =p cj2- 

We employ a numerical simulation in order to illustrate the dynamics of the 
system, dependent on the combination of laser parameters, as well as propagation 
effects of the coupled driving (fundamental), weak convertible and generated elec- 
tromagnetic waves in an optically-dense medium. The spacial behaviors of both 
sum-frequency and difference-frequency processes are considered, while the funda- 
mental radiation and consequently medium properties may vary along its length. 

2 Principal equations 

A partial energy-level scheme relevant to the above described four-wave-mixing 
process Wzp = 2loi ^u>2 controlled by a strong auxiliary infrared off-resonant field Est 
is depicted in Fig. 1. Levels m and I give the major contribution to the corresponding 
two-photon quantum pathways, and level / to Stark shift. When analyzing sum- 
frequency process at lu + , level I is placed above n with the same three-photon 
detuning. All the fields are assumed pulsed, where the duration of each pulse is 
much shorter than the shortest relaxation time in the quantum system. Through 
a time-dependent laser- induced Stark shift of level n, i.e., SCRAP, the field Est 
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Figure 1: Energy levels and coupled fields. 



allows the control of two-photon excitation by the field E\. By this process, the 
entire population of the ground level can be moved to the upper state n by the end 
of the pulse E\. Alternatively, maximum coherence at the two-photon transition 
gn can be achieved by making the populations of states g and n equal by the end of 
the pulse E±. Furthermore, with a judicious delay of the pulse E2, one can achieve 
maximum nonlinear polarization at the frequency uj^ and therefore high conversion 
efficiency of the visible radiation to the VUV range. 

2.1 Equations for interacting electromagnetic fields 

The propagation of optical waves in a non-conducting medium is described by the 
wave equation 

V 2 £-i^ = ^ (1) 
c 2 <9i 2 c 2 dt 2 ' y ' 

where V is the electrical polarization, and we assume the nonlinear medium to be 
isotropic and all the waves to be identically polarized and to propagate along the 
z-axis: 

£j(z, t) = Ke{Ej(z, t) exp[i(ujjt — kiz)]}, 

Vj{z,t) = Re{P j (z,t)exp[i(ui j t - k[z)]}. (2) 

Here kj is the modulus of a wave vector at frequency u>j, k'j is the wave vector of 
nonlinear polarization at the same frequency, and Ej(z,t) and Pj(z,t) are slowly- 
varying envelopes. 

Assuming a uniform field distribution over the cross section and taking into 
account Eqs. (2) and k w w/c, one can write Eq. (1) in the approximation of 
slowly- varying amplitudes as 
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where P = P L + p( FWM ) e ~ lAkz ^ an d Afc is the phase mismatch of the interacting 
waves. We restrict ourselves to the major process under consideration and assume 
phase matching to be fulfilled so that Ak = 0. Then by setting z' = z, t' = t — z/c 
and using 



dE dE dz' dE dt' _ dE IdE 

dz dz' dz dt' dz dz' c dt' ' 

dE dE dt' dE dz' _ dE 

~dt 'W~di + 'dz~'~dt ~~dt' 1 



we can express Eq. (3) in the form 

dE/dz' = 2nikP. (4) 
The polarization P can be calculated with the aid of the density matrix as 

P = Nptjdji + c.c, 

where N is number density of the resonance atoms, dji are electro-dipole transition 
matrix elements, and pij are density matrix components oscillating at the corre- 
sponding frequencies. The problem then reduces to calculating the off-diagonal 
elements of the density matrix. 

2.2 Density matrix equations 

In the case where the pulse durations of all the interacting fields are much less than 
any relaxation time in the system, all the relaxation terms may be omitted, and the 
density matrix equations take the form 



Pran ^ [^mnPn ~i~ GmgPgn] 3 

Pin = —i [GlnPn + GlgPgn] , 

Pgl = i [PgnGnl + PgGgl] , 

Pgm — i \PgnGnm Gg m pg] 
Pgn i[GgmPmn CjnnPgm 

+ Ggipi n — GlnPgl — Gfnpgf] 

Pn = 2 111! [G nm p mn + GnlPln + G n fPf n ] , 

Pnf = iG n fp n , p g f = ip gn G n f , p g = 1 - p n . (5) 

Here ptj = dpij/dt (hereafter we omit the prime on z and t, assuming a trans- 



formation to a new coordinate system), Gij(t) = —dijEk(t)/2h are the coupling 
Rabi frequencies, and dij are electric dipole transition momenta. In the coupled 
equations (5), we assume that initially only the lower level g is populated and that 
the populations of levels m, I and / are negligibly small during the entire process. 
We also neglect photo-ionization. Further, we assume two-photon quasi-resonancc 
£l gn = (2u>i — ojgn) *C u>ij and that the conditions for one-photon dctunings 

flij = (ujk - ujij) > dj » t" 1 (6) 
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are fulfilled, where Tj is the minimum pulse length. 

We look for the solution of the coupled equations in the form 



Pij(t) = rij(t)exp{-iQijt} , 
Gij(t) = 9k(t) exp {-iflijt} , 
G mn {t) = agi(t) exp {-i(fl mn t} , 



where r(t) and g(t) are slowly-varying envelopes, and a = d mn /d gm . Then using 
(6) and |fy <C |f2y-ry| for all the off-diagonal elements except r gn , we can rewrite 
the equations (5) as 

Tmn — [Q'9lT'n 5l^*<?fi] /Slrani 

Tin = [g2T n + g*Lr gn ] /ftln, 

Tgl = -[r gn 92 + r gg-} /Sigh 

r gm = - [a*g{r gn + gir g ] /fl gm , 

r gn = -i[gir mn - agiVgm + g-r ln 

- g2r g i - g* S t r 9f - Sl gn r gn ] , 

f n = 2 Im [a*glr mn + g^n n + gstTfn] , 

r n f = —gstfn/^nf, r gf = -r gn 9St/Slgf , 

r g = 1 - r n . (7) 



Here Sl n f = uj S t - w„/, Sl g i = 2wi - w 2 - u g i, Sl g f = 2wi + vst - u g f, etc. After 
simple algebra, Eqs. (7) reduce to 



21m 



9i 



SI, 



Slln 
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9gn 



Sin 
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(8) 



where 



9gn 



ag\ . 525- ag\ gig- 
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2.2.1 Generalized two-level scheme 

In the more general case of multi-level system a sum over all intermediate states m 
and I must be taken. In view of that, we introduce the values 



1 772 \9St\ n I TP |2 

n = -2- — =7i^i, s = — — = f3 s \Est\ 



i \gst\ 2 

s i=(^- + 7^-) |5i| 2 =/3|n|=/3i|^| 2 , 

\ * "gin * 'mn / 
.92.9- w p |.92 1 2 fl I 77 |2 

r-i = -2— — = ry 2 E 2 E-, s 2 = — — = ^l^l , 

s _ = \E-l = ,3_\E^, (10) 

where 7^ and /3j are polarizabilities depending on specific atoms and contributing 
quantum transitions. In the vicinity of a two-photon resonance fl m n ~ -U 9m and 
^in ~ — Qgi- Then Eqs. (8) take the form: 

r n = Im[(r* + r 2 )r gn ] , 

-i( ri +r 2 )(r n -l/2), (11) 

where 

n st = s + Sl + s 2 + s_. (12) 

It is seen from Eqs. (8) that £lst is the laser-induced shift of the two-photon 
resonance. 

2.3 Evolution of Rabi frequencies along the medium 

We assume that one-photon Rabi frequencies of all the fields vary in time as 

9i {t) = <?i exp{-(f - At,) 2 /2t 2 }. (13) 

Here 2n is the width of the i-th pulse (at e _1 of the power maximum), and An is 
the time delay of the i-th pulse relative to t = 0. The pulse widths are assumed 
to be much shorter than the relaxation rates of the system. We also assume that 
the spectral width of the pulse E\ is much greater than the Doppler width of the 
transition gn, and are replaced Eq. (4) by the equations for the corresponding Rabi 
frequencies, 

dgi/dz = -iirh~ 1 ki\d gm \ 2 N(r gm + a*r mn ) 

= -iai(r gm + a*r mn ), 
dg 2 /dz = -ia 2 ri n , dg_/dz = -ia^r gU 
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where oiifi — are resonant frequency-integrated absorption indices at the correspond- 
ing transitions for broadband probe radiations, with all other fields turned off. Then 
by scaling the medium length to the absorption length aZ , which is a characteristic 
parameter for the given medium that can be measured independently, we obtain 

dgi/d£ = -iKi(r gm + a*r mn ), 
= -iK 2 rin, 

dg-/d£ = -ir gU (14) 

where 

Ki = (fci|d gm | 2 )/(fc-M 2 ) 

= {knj g i\f gm \)/(k-u gm \f g i\), 
K 2 = (k 2 \d ln \ 2 )/(k-\d gl \ 2 ) 

= {k2U g l\fln\)/{k-t0ln\fgl\), 

£ = a-z. (15) 

In the above consideration we ignored Doppler effects, which is possible, for 
example, in atomic beams. In the case of warm gas, the solution of the equations 
must be averaged over a Maxwell velocity distribution of atoms. However, Doppler 
effects do not contribute much if the equivalent spectral width of the pulses or/and 
all detunings are greater or comparable with the corresponding resonance Doppler 
HWHM. 



3 Numerical simulation of SCRAP and maximum 
coherence in a generalized two-level scheme 

For further numerical simulations we introduce the parameters 5 = Q gn Ti, S = SoTi, 
R = ?*ioti, T = t/ri and St = Ar/ri, which are the detuning and maximum 
dynamic shift of the resonance scaled to the spectral width of the pulse E\ (Eqs. 
(13)), number of Rabi oscillations during the pulse duration n, and time and pulse 
shifts scaled to the duration of this pulse. In this subsection we analyze the dynamics 
of population transfer and induced coherence, assuming 92 = and the dynamic self- 
shift of the resonance to be negligible (/3 = 0). Then the equations (11) describe 
laser-induced oscillations in a two-level system controlled by the auxiliary Stark 
field. The corresponding level scheme and radiative coupling are shown in Fig. 2. 



3.1 Oscillations induced by a rectangular pulse 

First, we simulate the well-known dynamics induced by pulses of rectangular shape 
and duration r at Q gn = Q$t ( 7r /2, 7r and 2tt pulses). Then the equations (11) take 
the form: 



f gn = -«>i(r„ - 1/2), r n = -\ri\ 2 (r n - 1/2). 



(16) 
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Figure 2: Generalized two- level scheme and radiative couplings. 



Equation (16) describes induced oscillations with the two-photon Rabi frequency 
\ri\. The coherence r gn reaches its maximum {r gn = 0) at r n = 1/2. Then r g = 1 — 
r n = 1/2 too, and r gn , which is product of the corresponding probability amplitudes, 
reaches its maximum \r gn \ = 1/2 as well. 

At Q, = tt gn — £lst = and t = 7r/2|ri |, half of the population of the ground state 
transfers to the excited state, and the amplitude of the polarization oscillations \r gn \ 
reaches its maximum. At t = 7r/|ri| the entire population of the ground state is 
transferred to the excited state. At t = 3ir/2\ri\ the coherence again reaches its 
maximum, and at t = 2n/\ri | the system returns to its initial state. Such oscillations 
are illustrated in Fig. 3. Deviation from resonance leads to increased frequency and 
decreased amplitude of the oscillations for the level populations (Fig. 3 (left)), and 
the instant when maximum coherence is reached will vary. However, the maximum 
value of coherence, which is 0.5, is not achievable if the population of the upper 
level does not reach 0.5. (Fig. 3 (right)). 

3.2 Dynamics of generalized two-level system driven by Gaus- 
sian pulses 

In further analysis we shall account for the differences which appear for the pulses 
of Gaussian shape (Eq. (13)). This leads to variation in time of the two-photon 
Rabi frequency and Stark shifts as: 



Hereafter, we will scale a static deviation from the two-photon resonance Q gni am- 
plitudes of the dynamic resonance shift so and Rabi frequency rio to the spectral 



nit) = r i0 exp [-(i - An) 2 /rf] , 
s(t) = s exp[-(t- AT St ) 2 /T* t ] , 
st(t) =s l0 exp[-(<- Anf/rf] . 



(17) 
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Time T Time T 



Figure 3: Dependence of the population r n and coherence \r gn \ on time T in 
a generalized two-level scheme with excitation by a rectangular pulse at various 
deviations from resonance Sl gn - Time is scaled to the arbitrary value r which is 
much less then any relaxation rate in the system. The coupling Rabi frequency 
R = r\T = 27r/5, and S = Q, gn T are scaled to t _1 . Solid line -5 = 0, dashed - 
5 = 1.259, dash-dotted -8 = 2.5. 



width t 1 , and delay of the Stark pulse At and time t to the duration of the 
excitation pulse n: 8 = Cl gn Ti, S = son, R = riori, St = Ar/n, T = t/n. 

3.2.1 Persistent maximum coherence created by Gaussian pulses at neg- 
ligible dynamic self-shift of the resonance 

In the case under consideration in this sub-subsection (/3 = 0), the principal con- 
tribution to Qst is determined by Est, such that the entire Stark shift depends on 
time as 

n S t(t) = s cxp[-(<- AT Si ) 2 /4] ■ (18) 

We choose the intensity and frequency detuning of the Stark field such that the 
Stark shift fist can compensate for the two-photon detuning £l gn . In order to 
provide a dynamic two-photon resonance, the requirement |so| > \Q gn \ must be 
fulfilled. Then the resonance condition is met at two instants of time: 



<i, 2 = Ar st =F T St ydn(s AVi)- (!9) 

By a proper choice of the field parameters, one can provide first passage of the 
two-photon resonance at such an instant that the value of the two-photon Rabi 
frequency ensures a 7r/2 process over the period of passage through the resonance. 
For example, if the first instant of the crossing of the resonance is t\ = and the 
second one is £2 = 2Ars t , then the requirement is: 



Q gn = s cxp(-ArJ t /Tj t ) 



(20) 
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Figure 4: Population r« and persistent maximum coherence \r gn \ controlled by 
the dynamic Stark shift at various static deviations from the two-photon-rcsonancc 
and moderate two-photon Rabi frequencies. The Stark shift produced by the two- 
photon excitation is ignored. The duration of the control Gaussian Stark-shift pulse 
is rst/n = 1.6. Solid plot — 6 = 0, S = 0, R = 0.886; dashed — 5 = 5, S = 7.4, 
R = 3.18, St = 1.8; dash-dotted — S = 20, S = 23, R = 3.48, St = 1.34. 



A numerical analysis of corresponding feasibilities of creation of persistent max- 
imum coherence is presented in Fig. 4. Figure 4 (solid) illustrates the maximum 
coherence generated by the resonance ir/2 pulse with the Stark field turned off. 
Figures 4 (dash and dash-dotted) show such opportunities realized with the aid of 
the Stark field at various deviations of excitation from the two-photon resonance. 
The Stark shift at the front of the pulse is adjusted to compensate for the static 
deviation at the instant when the excitation is around its maximum values, while 
the Rabi frequency is about the magnitude required to transfer half of the popu- 
lation of the ground state to the excited one during the period of passage through 
the resonance. Due to the appropriately selected pulse shift, the Rabi frequency is 
too small for the population to return to the ground state while passing through 
the resonance at the rear of the Stark pulse. The results of simulations depicted 
in Fig. 5 show that the variation of the given radiative parameters about 10% of 
their optimum values leads to comparable change of the population transfer. As 
an important outcome, it demonstrates substantially less change of the coherence 
amplitude and therefore attests to its robustness. This is due to the fact that the 
product of the probability amplitudes of the ground and excited states changes less 
than the square modulus of probability amplitude of the upper state. 
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Figure 5: Robustness of maximum coherence created with the Gaussian pulses of 
two-photon excitation and Stark shift of the resonance, while the dynamic self-shift 
of the resonance is neglected (f3 = 0), and rst/n = 1.6. (a) - S = 5, R = 3.18, 
St = 1.8, solid - S = 7.4, dashed - S = 6.7, dash-dotted -S = 8.1; (b) - S = 7.4, 
R = 3.18, St = 1.8, solid line -6 = 5, dashed - 6 = 4.5, dash-dotted - 5 = 5.5; (c) 
- S = 23, S = 20, St = 1.34, solid - R = 3.48, dashed - R = 3.85, dash-dotted - 
R = 3.15; (d) - R = 3.48, S = 23, S = 20, solid - St = 1.34, dashed - St = 1.2, 
dash-dotted - St = 1.5. 
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Figure 6: Rapid adiabatic passage assisted by Stark-chirp and two-photon co- 
herence, while the dynamic self-shift of the resonance is ignored (J3 = 0), and 
rst/n = 1.6. (a) - 5 = 24, S = 75, St = 1.7, solid line - R = 30, dashed - R = 15, 
dash-dotted - R = 60; (b) - S = 75, R = 30, 8 = 24, solid - St = 1.7, dashed 
- St = 1.1, dash-dotted - St = 2.3; (c) - S = 24, R = 30, St = 1.7, solid line - 
S = 75, dashed - S = 50, dash-dotted - S = 150; (d) - i? = 30, S = 75, St = 1.7, 
solid line - 5 = 24, dashed -(5=12, dash-dotted - 5 = 48. 
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Time T Time T 



Figure 7: Effect of dynamic self-shift of two-photon resonance S\ — /3R on creation 
of persistent maximum coherence. S — 0. R = 0.886, solid line - /3 = 0, S = 
(7r/2-pulse), dash-dotted - (3 = — 1, S = 0, dashed (overlaps the solid one) - (3 = 
— 1, 6 = —0.65. Insert - dependence of the dynamic self-shift of the two-photon 
resonance si on time. The static deviation from the resonance is marked by the 
straight line. 



3.2.2 Stark-chirped rapid adiabatic passage and transient maximum co- 
herence 

As indicated above, the key idea of SCRAP [25] is to enhance the transition proba- 
bility at t\ and to minimize it at ti (see Eq. (19)), which is possible if the evolution 
is adiabatic at t\ and diabatic at t-i. Creation of maximum coherence is ensured 
in the process of rapid adiabatic passage of the entire population from the ground 
to excited state, where maximum coherence is always induced at instants when 
populations of the level becomes equal. Figure 6 presents such behavior. It shows 
robustness of such passage, although the instants when coherence reaches its max- 
imum may change substantially. 

3.2.3 Maximum coherence and rapid adiabatic passage controlled by 
the dynamic shift of a two-photon resonance produced by both 
fundamental and Stark Gaussian pulses 

In a multi-level system with appreciable difference between the principal quantum 
numbers of the ground and excited states, two-photon excitation may substantially 
contribute to the dynamic shift of a resonance. This makes transient effects and 
consequently the choice of the parameters ensuring maximum coherence and rapid 
adiabatic passage more complicated. Figure 7 presents such behavior, where the 
solid line shows feasibility to generate persistent maximum coherence in the absence 
of dynamic shift of the resonance (a 7r/2 pulse). It is seen that even a relatively small 
self-induced dynamic shift of the two-photon resonance may substantially change 
the time-behavior of the system (dash-dot plot). However, persistent maximum 
coherence can be restored through an appropriately chosen static deviation from 
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Figure 8: Population r n and coherence \r gn \ determined by the dynamic shift of 
the two-photon resonance driven both by resonant excitation s\ = (3R and by an 
auxiliary Stark-shift field s. R = 20, S = 0, f3 = -0.5; (a)- 6 = 0, (b)- 6 = -3, (c) 
- dashed line - 6 = —10, solid - S = — 15. Insert - dependence of the dynamic shift 
Si induced by two-photon excitation (dashed line) and by two-photon excitation in 
cooperation with the auxiliary Stark field s\ + s (solid line). 
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the resonance provided that the Rabi frequency is around the value corresponding 
to a 7r/2-pulse. 

Under larger intensities of two-photon excitation, the dynamics of the system 
controlled by the dynamic self-shift becomes more complicated, so that generation 
of persistent coherence without the means of independent control of the shift with an 
auxiliary pulse may become impossible (Fig. 8). Indeed, independent control of the 

0.8 



0.6 
0.4 
0.2 

-5 5 -5 5 

Time T Time T 

Figure 9: Creation of persistent maximum coherence controlled by appropriate 
static deviation from resonance and by appropriate delay of he auxiliary Stark- 
shift pulse in the system experiencing the dynamic self-shift of the two-photon 
resonance, si = (3R, (3 = —0.5, S = 15, S = 13.65, 6t = 1.6; dashed line - 
R = 15, solid - R = 30. Insert - dependence of the dynamic shift s% induced by 
two-photon excitation (dashed line) and by two-photon excitation in cooperation 
with the auxiliary Stark field si + s (solid line). The static deviation from the 
resonance is marked by the straight line. 

dynamic shift of the two-photon resonance with the aid of the auxiliary Stark-shift 
pulse, and based on the above discussed algorithm, ensures creation of persistent 
maximum two-photon coherence (Fig. 9). 

Figure 10 illustrates similar opportunities for transient maximum coherence gen- 
erated through persistent transfer of the entire population of ground state to the 
upper state (SCRAP). 

4 Dynamics of energy level populations, four- wave 
mixing at maximum coherence, and VUV gen- 
eration in Hg vapors 

In this section, we apply our theory for simulating the features of two-photon res- 
onant FWM controlled by SCRAP with the aid of a numerical model appropriate 
for quantum transitions in Hg. Then the energy levels depicted in (Fig. 1) are at- 
tributed to those of Hg as follows: g - 6s 2 1 5o, n - 8s 1 5o, m - 6p 1 Pi, I - 6p 1 Pi, 
f - 7p 1 Pi, and the laser parameters are similar to those described in Ref. [25]: 
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Figure 10: Persistent Stark-chirp-assisted rapid adiabatic passage controlled by the 
appropriately delayed pulses and by static deviation from resonance in the system 
experiencing a dynamic self-shift of the two-photon resonance, si = /3R, (3 = —0.5, 
S = 75, R = 20, 5t = 1.4; dashed line - S = 10, solid - S = 24. Insert - dependence 
of the dynamic shift si induced by two-photon excitation (dashed line) and by two- 
photon excitation in cooperation with the auxiliary Stark field si + s (solid), the 
static deviation from the resonance is marked by the straight line. 



Ai = 268.8 nm, X st = 1064 nm, A = 532 nm, n = 3 ns. Then A_ = 179.8 nm, 
A+ = 107.3 nm, 6 nf ~ 5 = 8.9 X 10 3 , S nl = -2.2 x 10 4 and 5 gm = -2.4 x 10 5 (all 
detunings are scaled to t^ 1 ). The parameters a and Kj are taken as a — 0.345, 
K_ = 1,K 2 = 0.04, Ki = 0.67. We have assumed that other levels, which are 
not depicted here, do not contribute substantially to the coupling. The other pulse 
durations are chosen as: t 2 = Ti, Tgt/ri = 1.6. The delay between all these pulses 
can be varied. We introduce the parameters Gi = \gi\Ti/2n and Z = £ti/27t, which 
denote reduced by 2ir the number of Rabi oscillations during the pulse E\ and the 
medium length scaled to the resonant absorption length at the transition gl for 
the nonmonochromatic radiation with the spectral width Tj~ . The results of the 
simulations are presented below. 

As outlined above, in order to achieve a transfer efficiency close to unity, one 
must fulfill the adiabatic condition described in Ref. [25]). Through a proper choice 
of Goi , Gost , static detuning 5 and pulse delay Sr s t , we can ensure various values and 
evolution in time of the populations and coherence at the entrance to the medium. 
Figures 11-14 demonstrate the possible dynamics leading either to maximum pop- 
ulation transfer or to equal populations and maximum coherence. These pictures 
(along with Fig. 6) illustrate the robustness of these processes, regardless of the 
significant change in some of the parameters. 

The generated radiation depends strongly on the dynamics of the coherence and 
on the delay at which the pulse G 2 to be converted is applied. For the cases depicted 
in Fig. 11, one can choose either a relatively short interval, when populations of the 
levels n and g becomes equal and the coherence amplitude reaches its maximum, or 
a 'plateau' with lower but almost constant amplitude of the two-photon coherence 
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Figure 11: Dependence of the population r n and coherence \r gn \ on time at the 
entrance to the medium. Goi = 910, Go s t = 325. Solid line - 5 = 0.5, Sr s t = —3, 
dashed - S = 5.6, 5r s t = 2. Other parameters are given above. 
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Figure 12: Dependence of the population r n and coherence \r gn \ on time at the 
entrance to the medium. Go st = 325, S = 0.5, Sr s t = —3. Solid line - Goi = 910, 
dashed - Goi = 1173, dash-dotted - Gqi = 525. 



\r gn \ in time. Furthermore, in our simulations, if it is not indicated otherwise, we 
assume that Gost = 325 and at the entrance of the Hg cell, Goi = 910 (which 
corresponds to R = 15, S = 75), and the amplitude of the pulse G2 is chosen to 
be non-perturbatively small (G20 = 1-6 x 10~ 2 ). The delay of the pulse G2 is set 
as ST2 = 5, so that this pulse overlaps the plateau in the time dependence of the 
coherence \r gn \ at the entrance to the medium (Fig. 11, solid line), or St2 = 0, so 
that the maximum of the convertible pulse coincides with the maximum of the time 
dependence of the coherence (Fig. 11, dashed line). Figures 15(a-d) computed for 
difference- frequency processes show that the shape of the generated and convertible 
pulses may substantially vary along the medium, subject to the coherence dynamics 
and the time delay St at the entrance of the medium. It is seen that if the pulse 
G2 overlaps the plateau of the time dependence of |r ffra |, there is no significant 
transformation of the pulse shape along the medium (the corresponding plots in 
Fig. 15(c) coincide). On the contrary, if the maximum coherence is reached in a 
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Figure 13: Dependence of the population r n and coherence \r gn \ on time at the 
entrance to the medium. Goi = 910, 5 = 0.5, ST st = —3. Solid line - Gost = 461, 
dashed - Go st = 266. 
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Figure 14: Dependence of the population r n and coherence \r gn \ on time at the 
entrance to the medium. Goi = 910, Go s t = 325, <5 = 0.5. Solid line - <5r sf = —4, 
dashed - <5r sf = —3, dash-dotted - Sr st = —2. 



transient regime, and the maximum of the pulse Gi coincides with the instant when 
the populations of the levels r g and r„ become equal at the entrance to the medium, 
the change of pulse shapes is most significant (Fig. 15(b,d)). In the first case, all the 
parts of the pulse G2 are converted homogeneously. In the second case, the wings 
of the pulse G2 are not converted and enhanced since the shape of the pulse \r gn \ 
is substantially narrower than that of G2 . This leads to a sharpening of both the 
generated and convertible (i.e., amplified) pulses. The shape of the excitation pulse 
may substantially change along the medium too, although its energy has practically 
no change (see Fig. 16). 

Figure 17(a-c) shows the dependencies of the pulse energy 

(0 J IG^l 2 dr - J \G 0i (r)\ 2 dr m 2 \d ln \ 2 , . 

ph f\G 02 (r)\^dr ' o;.!^! 2 ' (AL> 



20 



S. A. Myslivets et al. 




Figure 15: Shapes of the generated IG-/G02I 2 (a,b) and amplified IG2/G02I 2 (c,d) 
pulses at Z = 3 x 10 6 . (a) and (c) correspond to the plots in Fig. 11 (solid line), 
with St2 = 5, and (b) and (d) to the plots at Fig. 11 (dashed line) with 5t2 = 0. 
The dashed line in plot (b) is the input pulse scaled to the maximum of the output 
one. For the plot (c). the shapes of the input and output pulses are similar. 
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Figure 17: Evolution the pulse energy e p u of the convertible and generated fields, 
scaled to the initial number of photons at convertible frequency u>2 along the 
medium, (a-c). (a) - down-conversion lu- = 2lu\ — u>2 and (b) - up-conversion 
lu + = 2u>i + ll>2, both correspond to the plots at Fig. 11 (solid line) whith 8t2 = 5, 
and (c,d) to down conversion, which correspond to the plots at Fig. 11 (dashed line) 
whith St 2 =0. (c) - energy e p h, (d) - peak pulse power W^ x . For up-conversion 
(b) plots for generated (solid line) and for convertible (dashed) fields are different, 
for down-conversion (a,c,d) their behavior is similar. 
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Figure 18: Evolution of the pulse shape W p h(T) of the generated and convertible 
radiations (coincide), with all parameters the same as in Fig. 17(c). 




Time T Time T 

Figure 19: Evolution in time of the population r n and coherence \r gn | at the medium 
length Z =lx 10 6 . Goi = 910, G Qst = 325. Solid line - S = 0.5, Sr st = -3, dashed 
-5 = 5.6, 5r st = 2. 

and Figure 17(d) the pulse power at the instances when it reaches its maximum 

W (i) = l<?i(T) I mas - IgojMJ^ ^2 Mini 2 , . 

\Go2(rW max co^l 2 ' 1 1 

along the medium. Both are scaled to the corresponding value for the convertible 
field E2 at the entrance of the medium. The plots (a,b) are computed for the case 
when, at the entrance to the media, the pulse E2 overlaps the plateau in the time 
dependence of the coherence \r gn \. The simulations show that in this case there is 
no difference between the plots for energy and power. On the contrary, in the case 
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when the second pulse overlaps maximum of the transient coherence, the difference 
between energy and power becomes significant (Fig. 17(c,d)), because pulse shapes 
experience change along the medium (Fig. 18). 

The figures show that the evolution of the fields along the medium is strongly 
determined by the coherent nonlinear coupling. The difference in evolution of the 
field E2 in plots (a,b) is due to the fact that the field E% is enhanced in the difference- 
frequency mixing process and, by contrast, depletes during the course of the sum- 
mixing process. The simulation reveals a substantial difference in the evolution, 
depending on such coupling parameters at the entrance to the medium, like the 
static detuning and delays between pulses, which control the coherence in time 
and along the medium. The maximum achieved amplitude of the generated field 
may change by several times, and the conversion efficiency for the sum-frequency 
process may approach unity. The figures prove the feasibilities of substantial (by 
more than several tens times) improvement of the conversion efficiency by a judicious 
fit of the parameters at the medium entrance. Since the pulse shape of the driving 
field E\ varies along the medium, consequently, the populations r g , r„, as well 
as the coherence r gn , vary along the medium too. Therefore, the system being 
optimally prepared at the entrance (Fig. 11) may evolve to less optimum along the 
medium (Fig. 19, solid line). Alternatively, this can be properly adjusted so that the 
generated radiation may grow over a greater length, and the conversion efficiency 
may become substantially larger as well. Consequently, conversion efficiency of the 
input green radiation to VUV range, which is greater than unity, becomes possible 
for difference-frequency process. 

Simulations show that, besides other parameters, the actual power, energy and 
photon conversion efficiencies are strongly dependent on the strengths of the in- 
volved transitions and on their ratios (parameters a and Kj). Figure 20 displays 
the same results as in Fig. 17(a) except that the ratio \di n /d g i\ = 0.35 is changed for 
1.1. Consequently, the parameter K2 changes from 0.04 to 0.4. The squared mod- 
ulus of the nonlinear polarization, which determines the FWM generated power, is 
proportional to the product of squared modules of these transition electric dipole 
matrix elements. Consequently, the output generated power increases as well. (The 
squared modulus of the output Rabi frequency, plotted in some of the graphs, in- 
creases by several orders even at the same output power.) Therefore, the major 
features of the generation, optimum conditions and and maximum output power 
depend on the specific quantum system employed. 

The characteristic absorption length zq = (a_Ti/27r) -1 for Hg and the medium 
lengths corresponding to Figs. 17-20 are estimated as follows: The frequency- 
integrated absorption index is calculated as a_ = irT g iaQ^ = X 2 s giAi g N / '4, where 
«o- is the resonant value of the absorption index for the monochromatic radiation, 
r g ; is the half-width of the resonance, Ai g is the spontaneous relaxation rate at 
the transition, gi is the degeneration factor for level I, and N is the Hg number 
density. The oscillator strength for the Hg transitions 6s 2 1 So — &p 1 P\ (A = 184.9 
nm) is / ~ 0.96. Since Ai g = 0.67g g fi g /giX 2 g , where A is given in cm, one obtains 
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Figure 20: Same as in Fig.l7(a), but \di n /d g i \ = 1.1 (K 2 = 0.4). 



a- = (0. 67 g g fi g N/ 4) cm -1 c _1 , and for t = 3 ns, and N = 10 16 cm~ 3 , an estimate 
gives a s Ti/27r « 7.6 x 10 5 cm -1 . Therefore, Z = 10 6 corresponds to a medium 
length of about 1.3 cm. 



5 Conclusions 

We have presented an investigation of specific features of four-wave mixing of de- 
layed short pulses in two-photon resonant media, controlled by a dynamic Stark 
shift of the resonance. Through numerical simulation, we have analyzed different 
regimes of preparation of laser-induced maximum coherence and its implementation 
for frequency conversion. We have shown that by a proper choice of the time delays 
between the interacting laser pulses, the conversion efficiencies may be enhanced 
by several orders of magnitude with respect to conventional four-wave mixing. The 
results indicate the feasibility of this technique to generate strong tunable vacuum- 
ultraviolet radiation. 
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